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A Relationship Between off-Specular Scattering and 
Time Spread Phenomena 


1. Introduction 


Acoustic signals received from near sea surface explosive 
detonations were observed to spread in time of arrival due to 
bottom interaction. The time spread extends considerably beyond 
the pulse width. In this study, off specular scattering due to 
bottom roughness which is thought to be a possible mechanism 
causing the spreading, is investigated. The scattering theory of 
Marsh-Kuo -1 is utilized in an attempt to predict the observed time 
spreading. 

2. Selected Geometry 


For many experiments the source and receiver are 
approximately at the same distance K from the bottom. For this 
geometry, off-specular arrivals at the same time, t , are 
reflected from an elliptical ring: 

<&/!>$■" 1 , ( 

where ( 


■=. specular arrival time, 
ia - direct arrival time. 


C = acoustic speed. 


3. Model Description 

Utilizing a method of perturbation, Kuo* generalized the 
scattering theory of Marsli and predicted reflected field 
characteristics from a rough transmitting bottom. To the second 
order of perturbation, the reflected field consists of specular 
and scattered components. 

The specular component is a concentrated beam in a forward 
specular direction. Its energy level depends on Rayleigh's 
reflectivity and the amount of scattered energy lost in all other 


3 



TM 861175 

directions. The scattered energy lost in all directions depends 
on the entire range of the bottom roughness spectrum. 

The scattered component is diffused and its energy level in 
any direction is dependent on a particular part of the bottom 
roughness spectrum. Therefore, the scattering process in an 
off-specular direction is roughness wave number selective. 

Given locations of the point source and receiver as indicated 
in Figure 1, each ray emerging from the source is assumed to be 
approximately plane at the scattering boundary. The first 
arrival after the direct arrival is the specular component which 
is reflected at the boundary point with coordinates (0,0,0). The 
succeeding arrivals, at a particular time, are all off-specular 
components arriving simultaneously from a particular elliptical 
ring depicted in the same figure. Figure 1 indicates one such 
off-specular component arrival. To find the received level at a 
particular time, one must sum up all contributions from such 
arrivals along a particular elliptical ring. 

Kuo's^ scattering strength, m", predicts off-specular 
components to be proportional to the product of the scattering 
function and the roughness spectral level at a selected wave 
number. The scattering function is independent of the roughness 
but is dependent on the incident/reflected wave and media 
properties. The selected wave number is the difference wave 
number between incident wave number t$, =. (zir (i ) 

and the reflected wave number = C2iri/c ) < \,>*0 
where ^ = acoustic frequency, C*,(3> = incident wave 

direction cosines with respect to ( x, ^ ), ( )= reflected 

wave direction cosines with respect to ( ). The magnitude 

of the selected wave number can be expressed in terms of 
geometrical quantities as given by the following expression. 


tirM® i, [ £ viV*. (rf “ 

(4) 

+ * u***' T- 

where y-i. yf 

(5) 


(6) 

- ”217^ /C 


x, y: satisfies equation (1). 



The selected wave number given by (4) will vary from zero at 
specular point (0,0,0) to 2k^ at eo . The latter selected wave 
number is also the selected wave number in the back-scattering 
direction. The extent of variation in the selected wave number 
along each ellipse was investigated. It was found that the 
magnitudes of the selected wave numbers along any ellipse vary as 
much as a factor of 2 within the area between the source and 
receiver. However, outside of that area variations were small. 
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Received level (RL) can be expressed in terms of the 
scattering strength level, 10 log m", by the following 
expression: 


RL * 10 log m" + SL - TL + 10 log A 


(7) 


,where SL = source level, 

TL = total transmission loss level, 
A = insonified area. 


For negligible water medium absorption loss, the total 
transmission loss becomes the double spherical spreading loss 
given by the following expression. 


TL = 20 log r^ + 20 log rj 


( 8 ) 


The received level at a given time, , is considered to be 
contributed to by a sonified elliptical ring area of one pulse 
width. The corresponding selected wave number contained in the 
scattering strength varies over this insonified ring area and 
the effect of each elementary area on RL should be integrated 
over the entire elliptical ring. However, for this study, a much 
simpler approximate method described below is utilized. 

4. Results for a Rigid Boundary 

The mathematical model presented by Kuo^ in 1964, was 
intended for scattering from a non-shear supporting sedimentary 
layer of an infinite thickness. As such, the model can not 
really predict scattering from basement material. A model 
development for shear supporting sediment and basement layers is 
underway. One approximation that can be made at this time is to 
let the density of the sediment become very large in comparison 
with that of water. Effectively this causes the sediment to 
approach a rigid boundary and simulates solid basement material 
to some extent. 

For this approximation, the scattering strength (m") is 
simplified to the following expression: 


1 * T V ) U,-*r| 
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where E represents a boundary roughness spectrum. In deriving 
the above expression, bottom roughness is assumed locally 
isotropic n^ar the selected wave number region - on the order of 
0.1 (meter ). Spofford/Holmes/Spiess J show measurements 
indicating approximate isotropic roughness down to wave number of 
about 10 (meter A ). For the following calculation, the power 
law spectrum, 

E(k) = E Q k ,m = 0 to 3, will be assumed. 

Because of geometrical symmetry and symmetry of , and 
in equation (9), scattering strength varies only in one quadrant. 
Scattering strengths at three locations in one quadrant are 
calculated and averaged. The averaged value is then used in 
equation (7) as applicable to the entire elliptical ring area. 
Source level is assumed zero dB. 

Figure 2 depicts the time spread in received levels 
(propagation loss) given by off-specular scattering. Signal 
frequency and pulse length are 500 B and 20 ms, respectively. 

The specular arrival time, is calculated to be 13.849 seconds. 

The first off-specular arrival level is calculated at 30 ms after 
the specular arrival. The succeeding off-specular levels are 
calculated at 20 ms intervals. Each curve represents the result 
of using a different power law for the roughness spectrum. 

In the same figure, the received specular path loss, based 
only on the spherical spreading, is 86 dB. 

Figure 2 predicts a higher received off-specular level and a 
faster decay rate for a larger spectral power law. Over the 
range of arrival time, the selected wave numbers are in the range 
of about 0.09 to 0.3 (meter x ), i.e. magnitudes of less than one. 
The selected spectral shape E(k) = E k m , indicates that the 
case of higher power law assumes a h?gher spectral level due to 
the magnitude of k being less than one. Therefore, the greater 
received level for a larger m value is due to the higher spectral 
level used. For the same spectral level, E should be adjusted 
for different m. However, the higher decayrate for larger m is 
physically caused by the faster drop in the roughness spectral 
level as the arrival time increases which increases the magnitude 
of the selected wave number. 

3 

According to measurements of Spofford et al , the value of m 
is about 3. For this power law, the off-specular received level 
is found to be independent of frequency. For m = 2, frequency 
dependence is depicted in Figure 3. 

5. Results for Sedimentary Boundary 

Time spread catised by off-specular scattering from a 
transmitting rough boundary is investigated here. ^For this case, 
the full scattering strength equation given by Kuo -1 is utilized. 
The general equation is such that an average over two quadrants 
is needed. Otherwise computations are similar to the rigid 
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boundary case. 

Figure 4 presents received levels from various kinds of 
bottom boundaries represented by various combinations of bottom 
to water density ratio / and wave number ratio k. A 
sediment represented by ( f, -fe, ) = (1.7, 0.8) was previously 
used by Kuo because it presented the best fit to earlier 
experimental back-scattering data. This sediment type together 
with a basalt type bottom, gives received off-specular levels 
which are very close to those of the rigid boundary investigated 
previously. A sediment Soup type bottom (expected in the area 
far away from continental shelf) gives very low off-specular 
levels. 

Figure 5 presents frequency dependence. Though time decay 
rate is higher in general, frequency dependence is similar to 
that of the rigid boundary case. 

Figure 6 depicts dependence on the grazing angle. In all 
cases, the total transmission loss is geometrically adjusted to 
have the same value. 

6. Concluding Remarks 

The results presented here represent a first attempt to model 
the time spreading due to scattering using perturbation theory. 
Detailed comparison with experimental data is not performed. 
However, private communication with experimentalists do indicate 
that the predicted decay rate is representative for: some cases. 
There are some experimental data which show a much slower decay 
rate. For these cases further studies are needed. 

The possible causes of a much slower decay rate are many. 

A few possibilities are suggested below. 

1) Multi-rough boundaries associated with multi-layers, can 
cause multi-off-specular returns and cause a slower decay rate 
compared to a single rough boundary. 

2) The combination of boundary and volumetric scattering can 
also cause a slower decay. 

3) A bottom material that supports shear can alter both the 
decay rate and the off-specular scattering levels considerably. 
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